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CHARACTERIZATION OF THE CUBIC EXPONENTIAL FAMILIES 
BY ORTHOGONALITY OF POLYNOMIALS 

By Abdelhamid Hassairi and Mohammed Zarai 

Sfax University 

This paper introduces a notion of 2-orthogonality for a sequence 
of polynomials to give extended versions of the Meixner and Feinsilver 
characterization results based on orthogonal polynomials. These new 
versions subsume the Letac-Mora characterization of the real natural 
exponential families having cubic variance function. 

1. Introduction. Let F = {P(m,F),m £ Mp} be a natural exponential 
family (NEF) on the real line parameterized by its domain of the means 
Mp. If Vp{m) denotes the variance of the probability distribution P(m,F), 
then the mapping m Vjr(m) is called the variance function of the fam- 
ily F. The importance of the variance function stems from the fact that it 
characterizes the family F within the class of all natural exponential fam- 
ilies. Furthermore, for many common NEFs the variance function takes a 
very simple form. Morris (1982) describes the class of real NEFs such that 
the variance function is a polynomial function of degree at most 2 in the 
mean. Up to affine transformations and powers of convolution, this class 
includes the normal, Poisson, binomial, negative binomial, gamma and a 
sixth family called hyperbolic cosine and nothing else. The Morris class of 
quadratic NEFs has received a deal of attention in the statistical literature 
and many interesting characteristic properties have been established. We will 
be concerned here only with the properties based on the notion of orthogo- 
nal polynomials. A remarkable characteristic result is due to Meixner (1934) 
[see also Letac (1992)]. It characterizes the distributions \x for which there 
exists a family of /U-orthogonal polynomials with an exponential generating 
function. These distributions generate exactly the Morris class of NEFs. A 
second characterization is due to Feinsilver (1986), who shows that a certain 
class of polynomials naturally associated to a NEF is /i-orthogonal if and 
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only if the family is in the Morris class. In the present paper, we will be 
concerned with the class of real natural exponential families having cubic 
variance function (i.e., a variance function which is a polynomial of degree 
less than or equal to 3). In fact, Letac and Mora (1990) have extended 
the work of Morris by classifying all real cubic natural exponential fami- 
lies. They have added to the Morris class six other types of NEFs which 
may be obtained from the Morris class by the action of the linear group 
GL(R 2 ) [see Hassairi (1992)]. The best known among such families is the 
Inverse-Gaussian family. Our aim is then to extend to the Letac-Mora class 
of cubic NEFs the different characterizations established for the Morris class 
and based on orthogonal polynomials. We show in fact that the cubicity of 
the variance function is characterized by a property of orthogonality which 
will be called the property of 2-orthogonality. The result is interesting in its 
own right and seems potentially very useful for asymptotic expansions and 
approximation. In Section 2, after a review of exponential family theory, we 
specify some facts about the Feinsilver sequence of polynomials associated 
to a NEF; in particular, we show that the generating function of this se- 
quence converges in a neighborhood of 0. In Section 3, we state and prove 
our main result concerning the characterization of the Feinsilver sequence 
of polynomials corresponding to a distribution generating a cubic natural 
exponential family by a property similar to orthogonality. We also show 
that this sequence is characterized by a four-term recurrence relation, while, 
as it is well known, a sequence of orthogonal polynomials satisfies a three- 
term recurrence relation. In Section 4, we first determine the families of 
2-orthogonal polynomials with exponential generating function. [The fami- 
lies with a formal exponential generating function are sometimes called in 
the literature Sheffer polynomials of type 0; see Sheffer (1939) and Rainville 
(I960).] This leads to another characterization of the Letac-Mora class of 
cubic NEFs which may be considered as the extension to this class of the 
Meixner characterization. 

2. Exponential families and associated polynomials. We need first to 
review some facts concerning natural exponential families and to introduce 
some notations. If fj, is a positive Radon measure on the real line R, we 
denote by 



its Laplace transform, and we denote by 0(/i) the interior of the convex set 
D([i) = {6 G R;L M (0) < oo}. A4(R) will denote the set of measures \x such 
that 0(/i) is not empty and /i is not concentrated on one point. If /i is in 
.M(R), we also denote 




k ll (6) = lo g L fM (e), 0e ©(//), 
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the cumulate function of fi. 

To each [i in .M(M) and 9 in 0(/i), we associate the following probability 
distribution on R: 

(2.1) P(0, n)(dx) = exp(9x - k^(9))fi(dx). 
The set 

F = F( f ,) = {P(e, fJ ,y,eeQ(ri} 

is called the NEF generated by /x. If \i and // are in M (R) , then F{y) = F(fi') 
if and only if there exists (a,b) in M 2 such that /i'(dx) = exp(ax + b)fi(dx). 
Therefore, if /i is in A4(R) and F = F(/j,), 

B F = {fi'eM(M.y,F(»')=F} 

is the set of basis of F. 

We have, for all 9 in B(/i), 

(2.2) = J R X n eMOx)Kdx) < +oo, 

and consequently, all the moments of P(9, jj) are finite. In fact, 

The function is strictly convex and real analytic. Its first derivative 
defines a diffeomorphism between 0(/i) and its image M F . Since k'^9) = 
J^xP(9, fj,)(dx), M F is called the domain of the means of F. The inverse 
function of k' is denoted by Vv an d> setting P(m,F) = P(ip(m), //), the 
probability of F with mean m, we have 

F = {P(m,F);me M F }, 

which is the parameterization of F by the mean. 
The density of P(m, F) with respect to /U is 

(2.4) f^(x,m) =exp{^ M (m)x-fc M (^ M (m))}. 
For m in Mp, we denote 

(2.5) V F (m) = ( (x-m) 2 P(m,F)(dx). 

JR 

Then 

(2.6) V F (m) = k'^(m)) = (^(m))" 1 , 

and, the map m \— > V F (m) is called the variance function of F. It entirely 
characterizes the NEF; that is, if F and F' are two NEF such that V F (m) = 
V F '(m) on a nonempty open set included in M F n M F >, then F = F' . 



(2.7) -^—^f fl (x,m) = P n (x,m)f f ,(x,m) 
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Consider now a real natural exponential family F and take [i = P(m ,F) 
with m fixed in Mf- The density / M (-,m) of P(m,F) with respect to fi is 
still given by (2.4) with /J^mJ = 1. It is easily verified by induction on n 
in N that there exists a polynomial P n in x of degree n such that 

gn 

dm n ' 
and 

(2.8) P n+ i(x,m) =^' ll (m){x - m)P n (x, m) + R n +i(x, m), 

where R n +i is a polynomial in x of degree < n + 1. In particular, we have 
that 

(2.9) P (x,m) = l and Pi(x, m) = ip' fM (m)(x — m). 

We now make a useful observation through the following theorem whose 
proof will be given in Section 5. For the sake of simplification, we set 

(2.10) P n (x) = P n (x,m ). 

Theorem 2.1. Let F be a NEF on R and let fi be a fixed probability in 
F with mean m . Let P n (x) be the polynomials defined by (2.7). Then 

E '^^^ 

is an entire series in L 2 (fi) of nonzero radius of convergence. 

It should be remarked that there exists r > such that, for all m in 
]m — r, m + r [ and for all x in R, 

(2.11) /,(x,m) = ^ (m ~7° )W p w (x). 

», fix 

To conclude this section, we mention that in many interesting situations, 
P n can be calculated by mean of the Faa di Bruno formula 

where k = k\ H h k n and the sum is taken for all integers kj > such that 

k\ + 2A;2 + • • • + nk n = n. That is, if we denote 

a{k){x) = ^? [( ^ (m)x " MVv(™)))L= mo > 

then 

[ '~ , nt ^ , h\---k n \[ 1! ) "A n! J 

fci+2fc 2 H \-nk n =n Tl \ / \ / 
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with kj in N, for 1 < j < n. 

The most famous example in this topic is the inverse Gaussian distribution 
with parameters 1/2 and p > defined by 

K dx ) = exp \ 2x~ J ho.+ool^ dx - 

The NEF generated by \i belongs to the Letac-Mora class. In fact, a standard 
computation shows that 6(u) =] — oo,0[ and k /1 (6) = —p\J—2Q. If F = F(/j,), 
then Mp =]0,+oo[, ipJm) = —p 2 /2m 2 and Vp(m) = m 3 /p 2 . 



For all m £ Mp, we have 



fci+2fc2H hifc n =n 



! ,„/ 2 1 \ kl /n + 1 1 



-7TX- 



3. Characterization of the cubic families in the Feinsilver way. As pointed 
out in the Introduction, any real cubic natural exponential family can be ob- 
tained from a quadratic family via the action of the linear group GL(M. 2 ). For 
instance, the quadratic natural exponential families on R have been charac- 
terized by Feinsilver (1986) as the ones for which the polynomials P n (x) are 
u-orthogonal. In this section we show that the polynomials P n (%) associated 
to a cubic natural exponential family have also a characterizing property of 
orthogonality which will be called the property of 2-orthogonality. 

Definition 3.1. Let u be a measure on R such that J \x\ n fi(dx) < oo for 
all n £ N. A family (Q^neN of polynomials on R is f/, — 2-orthogonal if, for all 
n and q in N*, / Q n (x)Q q (x)/j,(dx) = when n > 2q, and / Q n (x)/j,(dx) = 0. 

Next we give our first main result. 



Theorem 3.1. Let F be a NEF on R and let u be an element of F with 

d n f ( 

9m" 



mean m . Consider the polynomials (P n ) n gN defined by P n (x) = -jf-wfn(x, 



m)\ m=mo . Then the three following statements are equivalent: 

(i) The polynomials (P n ) ng pj are fi — 2-orthogonal. 

(ii) F is cubic. 

(hi) There exist real numbers (aj)o<i<3 such that, for all n > 2, 
xP n {x) = a 3 A 3 n P n _ 2 {x) + n(a 2 {n - 1) + l)P n _i(x) 
+ (nai + m )P n (x) + a P n+ i(x), 
with = n(n — l)(n — 2). Furthermore, in this case we have 

Vir(m) = 03(111 — m G ) 3 + 02(171 — m ) 2 + a\(m — m ) + a Q . 



G 
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To help in the proof of this theorem, let us give in Table 1, for each 
of the six types of NEF on K with polynomial variance function of degree 
3 [see Letac and Mora (1990)], the sequence of P(m ,F) — 2-orthogonal 
polynomials P n (x) defined by its recurrence relation for m = 1. 

Proof of Theorem 3.1. (i) =>• (ii). Equation (2.11) allows us to claim 
the existence of r > such that, for all m in the interval }m — r,m + r[ 
and for all x in R, 

/,(x,m) = ^ (m ~7° )W p w (x). 

neN 

If, for (m,mf) E (]m — r,m + r[) 2 , we set 

(3.1) g(m,m') = exp{fe At (^ /Lt (m) +tp^(m')) - k^ip^m)) - k^tp^m'))} , 

then the /i — 2-orthogonality of the polynomials (P n ) and Theorem 2.1 imply 
that 

g(m,m') = J f^(x,m)f fl (x,m')n(dx) 

f sr (m-m ) q (m' -m ) n 

= / Ti P n {x)P q {x)n{dx) 

J n,qen n - q - 

x - (m — m ) q (m' — m ) n f _ . . _ , . , , . 

x (m — m ) q (m' — m ) n „ , , 

= i + X) ^ / ^»(*W*H«*«0. 

geN*,ne[(g+l)/2,2g-l] q ' 

Taking the derivative of (3.1) with respect to m, we get, for all (m,m') E 
(]m -r, m + r[) 2 , 

(3.2) = ^ ga„,(m-ra () )''" 1 (ra'-m () ) n , 

g>l,n6[(9+l)/2,2g-l] 

with a ng = ^/P n (x)Pg(x)^(dx). 

Making m = m in (3.2), then, since ijj^(m ) = 0, we get 

^(m )(m' - m ) = a n (m' - m ). 

This is true for all m! E ]m — r, m + r[; then 

a n = V^( m o)- 
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Table 1 



Type 



Induction relations 



Inverse Gaussian with parameter p = 1 
Vf(to) = (m — m ) 3 + 3m (m — m ) 2 
+ 3m 2 (m - 



+ m 3 



Strict arcsine with parameter p = 1 
Vf(?ti) = (m — m ) 3 + 3m„(m — m ) 2 

+ (3m 2 + 1) (m — m ) + ml + m c 



Takacs with parameters p = 1 and a ■ 
VF(m) = 2(m - m D ) 3 

+ (6m D + 3)(m — m ) 2 

+ (6m 2 + 6m + l)(m — m 
x (2m 3 +3m 2 +m D ) 



Large arcsine with parameters p = 1 
and a — 1 

^3 



V>(m) 



; 2(m — m D 
+ (6m D + 2)(m- 



m o y 



+ (6m D + 4m D + 1) (m — m a 
+ 6m 3 + 2m 2 + m D 



Ressel with parameter p = 1 
VHro) = 2(m - m ) 3 + (6m D + 3)(m - 
+ (6m 2 + 6m D + l)(m — m D ) 
x (2m 3 +3m 2 +m ) 



Abel with parameter p = 1 1 
Vjr(m) = (m - m D ) 3 + (3m D + 2)(m - m D ) 2 
+ (3m 2 + 4m D + l)(m — m Q ) 
x (m 3 + 2m 2 +m ) 



Po(x): 
P 1 (X)- 

P 2 (x)- 

P n +l(x)- 



Po(x): 
P 1 (X)- 

P 2 (x)- 

Pn+l(x) = 



P, 



P,(x) 

Pl(z) 

P 2 (x) 

(x) 



P„+i(x) 



1, 

(x-1), 



x z -6x + 3, 
(x-3n-l)P„(x) 
-n(3n-2)P„_ 1 (a;) 
-A 3 P„_ 2 (2;), ?i > 2 

1, 

i(ac - 1), 
i(a; 2 -5a; + 3), 

|[(ar-4n-l)P n (a5) 
-n(3n-2)P„_i(x), 



■i4»P„_ a (i)], 



n > 2 



= 1, 



Po(a: 
Px(x) = ±(x-1) 

P 2 (x) 
P«+i(x) 



15a; + ! 



|[(x-13n-l)P„(x) 
-n(9n-8)P„_i(x) 
-2A 3 P n _ 2 (a;)], 



n > 2 



= 1, 



Po(x) 

Pi(x) = ±(x-1) 
P 2 (x) 
P„+i(x) 



-L(x 2 - 13a; + 3), 



§[(x-lln-l)P„(x) 
-n(87i-7)P„-i(x) 
- 2A 3 n P n - 2 (x)], 



n > 2 



§(*-!), 
±(x 2 -7x + 4), 

i[(x-5n-l)P n (a;) 
-n(4ra-3)P„_i(x) 
-A 3 P n _ 2 (x)], 



n > 2 



Po(x) = l, 



Pi(x) = |(x-1), 
P 2 (x) = i(x 2 -15x + 8), 



|[(x-8n-l)P n (x) 
-n(5ra-4)P n _i(ar) 
- A a n P n . 2 (x)}, 



n > 2 



<s 
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Again we take the derivative of (3.2) with respect to m and we let m = m . 
We get, for all m! in ]m — r, m + r[, 

^'/i( m o)(m' - m ) + al 1 (V F (m') - V F (m )) + af^m' - m ) 2 
= 2a 22 (m' - m ) 2 + 2023(771,' - m Q ) 3 . 

Therefore, 

V F (m') = 2aj~ 1 2 a 23 (m / - m ) 3 + (2a 22 - ai 1 )a{ 1 2 (m' - m ) 2 

- a^ 2 tl)^{m )(m' - m ) + V F (m ). 

This implies that V F is cubic on ]m — r, m + r [ and, by extension, we obtain 
that F is a cubic NEF. 

(ii) =^ (hi). From (ii), there exist real numbers (aj)o<«<3 such that 

Vf = 0,3(171- m ) 3 + a 2 (m - m ) 2 + ai(m - m ) + a Q . 

On the other hand, we know that there exists r > such that, for all m G 
]m — r, m + r[ and for all x S R, 

Y j -2- P n( x ) = exp{V' M (77T,)x - k^ifj^m))}. 

nGN n ' 

Denoting 9 = ip^m), this may be written as 

(3.3) exp(fe) = ( YW ~ ^oT^f-) exp{fc M (0)}. 

V nGN n - / 

Taking the derivative with respect to 9 of (3.3) gives 
xexp(#x) = Y ^l( n (m - m ) n - l k''{9) + (m - m ) n k'(9)) exp{kJ9)}, 

nGN n - 

which is equivalent to 
Y ^xP n (x) 

nGN 

= ^ ^^(n(m - mo)™" 1 ^^) + (m - m D ) n m) 

nGN n ' 



E ^^f«5Za fe (m-m ) 

nGN \ k=0 



n+k-i + ( m _ mo ) ri + 1 + m {m - m ) T 



By identification, we get 

xP n (x) = a 3 n(n - l)(n - 2)P n _ 2 (x) 

+ n(a 2 (re - 1) + l)P n _i(x) + (nai + m )P n (x) + a P n+ i(x) 
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and (iii) is proved. 

(iii) =>■ (i). The result is easily obtained if we verify the three following 
facts: 

(a) For all n G N*, / P n {x)n{dx) = 0. 

(b) There exist real numbers /3* such that, for all n, g£N* verifying 
n>2q, 

x«P n (x)=(3l q P n _ 2q (x)+ Pn, q Ps(x), 

n—2q+l<s<n+q 

where (3^ q = if n = 2q \/q£ N*, 

(c) There exist real numbers (a q )o< q < n such that 



Pn(x) — Oi n X -\- ^ ^ QZqX^ . 

0<q<n-l 



Proof of (a). We first observe that 

/ ^•^( x ' m ) /i ( cte ) = V^O) f (x-rn)fn(x,rn)(i(dx) 

= ip'^frn) J (x — m)P(m, F)(dx) 
= 0. 

Since, for all n, we have 




1 



H(dx)= J \P n (x -m)\f IJL {x,m)ix{dx) 

= J \P n (x - m)\P(m,F)(dx) < +oo, 



[see (2.3) and (2.7)], then 



Ux,m)}fi(dx) = — / ^-zrfu{x,m)ii(dx). 



dm \ dm n J ^ ' J 9m 7 9m™ ^ 

r>m" 



Hence we obtain that, for all n G N*, / ■^^f fl (x,m)fi(dx) = 0. This, for 
m = m D , gives / P n (x)n(dx) = 0. 

Proof of (b). We can write (iii) as 
(3.4) xP n (x)=Pl l P n . 2 (x)+ Yl #.i p »( a 

n— l<s<n+l 

where i = a2,n(n — l)(n — 2). 
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For a fixed n in N*, let us show by induction that, for all q in N* such 
that 2q < n, we have 

(3.5) x0P n (x)=(3° ntq P n _ 2q (x)+ ]T P s n , q P.s(x), 

n— 2q+l<s<.n+q 

where @®,q = if n = 2g. 

For q = 1, it is nothing but equality (3.4). 

Suppose now that (3.5) is true for q and that 2(q + 1) < n. Then we have 
x q+1 P n {x) =x{x q P n {x)) 

= Pn,qXPn-2q(x) + 0n >g xP s (x) 
n—2q+l<s<n+q 

= 0nJ Pn,n-2qPn-2q-2(x) + £ ftUg.l^OO} 

I n-2q-l<s'<n-2q+l ) 



E f3 s n Jpl 1 P s -2(x)+ J2 rfaPA*)) 

-l<s<n+q I s-l<s"<s+l J 



+ 

n-2q+l 

= a 3 (3° q A 3 n _ 2q P n _ 2(q+1) (x) + Pn,q+l P s( x )- 

n-2(q+l)+l<s<n+q+l 

Hence there exist (^q+i) sucn that 

X q+1 P n (x) = (3l q+l P n . 2q (x) + Pn,qPs{x)> 

n—2q+l<s<n+q 

where 0° >ff+1 = a 3 [3^ q A 3 n __ 2q and [3^ q+1 = if n = 2(q + 1). 

Proof of (c). Since -Pi(x) = y F ^ m ) (ra — ?ti<o)> it is easy to show by in- 
duction that 

and this concludes the proof. □ 

4. Characterization of the cubic families in the Meixner way. This sec- 
tion is devoted to the characterization of the 2-orthogonal polynomials on 
M. with exponential generating function. 

We say that the generating function of the sequence of polynomials Qn 
is exponential if there exist r > and two real analytic functions a and b 
defined on ] — r, r[ such that, for all z in ] — r,r[, 

z n 

(4.1) y^Q n (x)— = exp{a{z)x + b(z)}. 
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Families (Q n /n\) satisfying (4.1) are considered in the literature under 
the name of Sheffer polynomials of type 0; see Sheffer (1939) and Rainville 
[(I960), Chapter 13]. Actually they are slightly more general since conver- 
gence in a neighborhood of is not required and (4.1) is considered as an 
identity between formal series in z. 

Theorem 4.1. Let F be a NEF on R and let y, be an element of F with 
mean m . Suppose that (Q n )ngN is a family of y — 2-orthogonal polynomials 
such that Q n is of degree n. Then the generating function of {Q n )nen is 
exponential if and only if there exists t S R* such that, for all n G N, 

Q n (x)=t n P n (x), 

where (P n ) is defined by (2.7). 

In this case, a(z) = ip^itz + m ) and b(z) = —k lJi {'^^{a{z)). 

Proof. Up to Q n = Q n /Qo, we can suppose Q Q = 1. 
<= Is obvious. 

=>■ There exist r > such that, for all z S ] — r, r [ , 



j[Y J Qn{x)-Ay{dx)= I \Y,Qn(?)Qo(x)—y 
J VneN n - ) J VneN 71 ' , 



n(dx) 



Q {x) 2 y(dx) 



= 1. 

On the other hand, writing the generating function of (Q n ) as in (4.1), we 
have 

/ C X] = J ex P{°( z ) x + K z )}K dx ) 

= exp{fe /J (a(z)) +b(z)}. 



.nGN 



Hence 

(4.2) b(z)=-k fl (a(z)). 
Proceeding similarly, we have that 

(4.3) j(^Q n {x)Q 1 {x)^ii(dx)=^j Q^xf^dx^z. 

Then Q\ is a polynomial of degree 1 in x. Therefore there exists u £ R* and 
v G R such that 

(4.4) Qi=ux + v. 
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Since / Qi(x)Q (x)fx(dx) = J Q\(x)fj,(dx) = 0, we get b= —um and 
Qi(x) 2 fi(dx) = J u 2 (x — m )(x — m )fi(dx) = u 2 Vp(m ). 
Furthermore, using (4.2)— (4.4), we get 

Qi(x) 2 fi(dx) ) z = / ex.p{a(z)x + b(z)}Qi(x)fi(dx) 



and we deduce that 



= u J (x — m ) exp{a(z) — kfj,(a(z))}fj,(dx) 

= U J (x — m )P(a(z),fj,)(dx) 
= u[k' (l {a{z)) -m ], 

u 2 VF(m )z = u[k' (a(z)) — m \. 



Therefore, k'(a(z)) = uVp(m )z + m , that is, a(z) = ip^(uVp(m )z + m ) 
and t = uVp(m ). 
Finally, we obtain 

V] Qn{x)—r = fjx, V F (m )uZ + TTl ). 

Corollary 4.1. Let F be a NEF on R and let fj, be an element of F 
with mean m . Then there exists a family of fj> — 2-orthogonal polynomials 
with an exponential generating function if and only if F is cubic. 



Proof. Follows easily from Theorems 3.1. and 4.1. □ 



5. Proof of Theorem 2.1. This section is devoted to the proof of Theo- 
rem 2.1. 

Let x be a fixed real number. The function m *—> fu(x,m) is real analytic 
in the interval Mp. Therefore, there exist an open set O of C containing Mp, 
an open set U of C containing @(/i) and analytic functions tpi and k\ such 
that ^>i|mf = t/V) &i|e(u) = hp an< i ^i(^) C U. For all x S R, the function 
z^ fi(x,z) = exp{xipi(z) — ki(ipi(z))} is analytic on ft. 

Since G @([J>) and m £ Mp, there exist a > and r > such that: 

(i) [-a, a] C 6(/i); 

(ii) the open disk D(m ,r) C f2; 

(iii) ViP(m o ,r))cZ>(0,f); 

(iv) V\ = (^i) _1 has no zero in D(m ,r). 
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Let us show that the function <j)(z) = fi(-,z) is well defined and continu- 
ously differentiable from D(m ,r) into L 2 (fj,). 
Since 



+00 



\fx(x,z)\ 2 »(dx) 



ex-p{x(ipi(z) + ) - (ki(ipi(z)) + ki(ipi(z) )}fi(dx) 



and 9 = ipx(z) + Mz) G [=fs %] C 9(/x), then 



3 ' 3 
+00 



\fx(x,z)\ fi(dx) < +00, 



and 4> is well defined. 

For the differentiability, if z Q is in D(m , r), then ( Jj/i)(-, z ) is an element 
of L 2 ([i) because 



x — z 



Vi(z) 



exp{x(ipi(z) + ipi(z) )}n(dx) < +00. 



We will verify that /l(,,Zo+ ^ h{ ' ,Zo) - j%fi(-,z) converges to in L 2 (n) when 
h converges to 0, that is, 

1 r + °° . , . , . , , . / d 



1 r 

(5.1) lim / 

fc->0 h 2 J-, 



fx(x,z + h) - fi{x,z ) -h[ — fx )(x,z 



fj,(dx) = 0. 



Writing the Taylor formula with integral remainder, 



f(h) - /(0) - hf'(0) = h 2 I (1 - u)f"(z + uh) du, 
(5.1) is equivalent to 







I r+00 rl 

(5.2) lim -rr I h 4 / —^fx(x,z + uh)(l-u)du 
h^O h A J-00 Jo oz z 



n(dx) = 0. 



But we have that 

d 2 



x-z \ 2 V x {z) + (x - z)V({z) 
Vi(z) 



V 2 (z) 



fi(x,z). 



Then Jp-/i(x, z Q + h)/f(x, z Q + h) is a second-degree polynomial in x whose 
coefficients are continuous in h. It is then bounded in all compact \h\ <h . 
Hence 



dz 2 



fx(x,z a + h) 



< \ax + bx + c\e ° 2 



with 



< 



2a 
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This with the dominated convergence theorem implies that 




r+oo rl Q2 2 

lim / / (1 — u)— —~fi(x, z + uh) du u(dx) 




+oo rl Q2 



\ C 1 ~ u )jr^h( x i z o)du n(dx)<+oo 



-oo JO 



and (5.2) follows. 

Hence <f> is continuously differentiable and so it is analytic on D(m ,r) 
and, in particular, if 
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